MATRIX EXTENSION WITH SYMMETRY AND ITS APPLICATION 

TO FILTER BANKS* 

BIN HANt AND XIAOSHENG ZHUANG^ 

Abstract. Let P be an r X s matrix of Laurent polynomials with symmetry such that P(z)P* (z) = 
Ir for all z G C\{0} and the symmetry of P is compatible. The matrix extension problem with 
symmetry is to find an s X s square matrix Pe of Laurent polynomials with symmetry such that 
[/r,0]Pe = P (that is, the submatrix of the first r rows of Pe is the given matrix P), Pe is parauni- 
tary satisfying Pe(z)Pe(^) = Is for all z G C\{0}, and the symmetry of Pe is compatible. Moreover, 
it is highly desirable in many applications that the support of the coefficient sequence of Pe can 
be controlled by that of P. In this paper, we completely solve the matrix extension problem with 
symmetry and provide a step-by-step algorithm to construct such a desired matrix Pe from a given 
matrix P. Furthermore, using a cascade structure, we obtain a complete representation of any r X s 
paraunitary matrix P having compatible symmetry, which in turn leads to an algorithm for deriving 
a desired matrix Pe from a given matrix P. Matrix extension plays an important role in many areas 
such as electronic engineering, system sciences, applied mathematics, and pure mathematics. As an 
application of our general results on matrix extension with symmetry, we obtain a satisfactory algo- 
rithm for constructing symmetric paraunitary filter banks and symmetric orthonormal multiwavelets 
by deriving high-pass filters with symmetry from any given low-pass filters with symmetry. Several 
examples are provided to illustrate the proposed algorithms and results in this paper. 
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1. Introduction and Main Results. The matrix extension problem plays a 
fundamental role in many areas such as electronic engineering, system sciences, math- 
ematics, and etc. To mention only a few references here on this topic, see [TJ [2l |3l O [BJ 
III [9j [121 El HSl [m [171 [19 . 20 . For example, matrix extension is an indispensable tool 
in the design of filter banks in electronic engineering ( [13l [HI [191 120] ) and in the con- 
struction of multiwavelets in wavelet analysis ([Il[l[3[l[Sl[Sl[71[i[Tni[Il[IIl[ISl[IS]). 
In order to state the matrix extension problem and our main results on this topic, let 
us introduce some notation and definitions first. 

Let p(z) — 'Ylik&PkZ^TZ G C\{0} be a Laurent polynomial with complex coeffi- 
cients pk G C for all fc G Z. We say that p has symmetry if its coefficient sequence 
{Pk\k£Z has symmetry; more precisely, there exist e G { — 1, 1} and c G Z such that 

Pc-k ^ spk, VfcGZ. (1.1) 

If e = 1, then p is symmetric about the point c/2; if e = —1, then p is antisymmetric 
about the point c/2. Symmetry of a Laurent polynomial can be conveniently expressed 
using a symmetry operator S defined by 

Sp{z):=-^ zGC\{0}. (1.2) 

When p is not identically zero, it is evident that (|l.ll) holds if and only ii Sp{z) — ez'^. 
For the zero polynomial, it is very natural that 50 can be assigned any symmetry 
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pattern; that is, for every occurrence of 50 appearing in an identity in this paper, 50 
is understood to take an appropriate choice of ez'^ for some e e 1} and c G Z so 
that the identity holds. If P is an r x s matrix of Laurent polynomials with symmetry, 
then we can apply the operator S to each entry of P, that is, SP is an r x s matrix 
such that [5P]j.fc :— S{[P]j^k)^ where [P]j.k denotes the (j, fc)-entry of the matrix P 
throughout the paper. 

For two matrices P and Q of Laurent polynomials with symmetry, even though 
all the entries in P and Q have symmetry, their sum P + Q, difference P — Q, or 
product PQ, if well defined, generally may not have symmetry any more. This is one 
of the difficulties for matrix extension with symmetry. In order for P ± Q or PQ to 
possess some symmetry, the symmetry patterns of P and Q should be compatible. For 
example, if SP = 5Q, that is, both P and Q have the same symmetry pattern, then 
indeed P ± Q has symmetry and S{P ± Q) = SP — SQ. In the following, we discuss 
the compatibility of symmetry patterns of matrices of Laurent polynomials. For an 
r X s matrix P{z) = X^fcez PkZ^ , throughout the paper we denote 

P*(z) := ^Pfe*z-'^ with P^:=pf, k € Z, (1.3) 

T 

where Pk denotes the transpose of the complex conjugate of the constant matrix Pk 
in C. We say that the symmetry of P is compatible or P has compatible symmetry, if 

SP{z) = {S9inz)S92{z), (1.4) 

for some 1 x r and 1 x s row vectors 9i and 02 of Laurent polynomials with symmetry. 
For an r X s matrix P and slu s x t matrix Q of Laurent polynomials, we say that 
(P,Q) has mutually compatible symmetry if 

SP{z) = {S9i)*{z)S9{z) and SQ{z) = {S9)* {z)Se2{z) (1.5) 

for some Ixr, Ixs, Ixt row vectors 9i, 9, 92 of Laurent polynomials with symmetry. 
If (P, Q) has mutually compatible symmetry as in (jl.Sp . then it is easy to verify that 
their product PQ has compatible symmetry and in fact 5(PQ) = {Sdi)*Sd2- 

For a matrix of Laurent polynomials, another important property is the support 
of its coefficient sequence. For P — J^kez, ^kZ^ such that Pk — ^ for all k £ Z\[m, n] 
with P„j 7^ and P„ 7^ 0, we define its coefficient support to be coeffsupp(P) :— [m, n] 
and the length of its coefficient support to be |coeffsupp(P)| := n — m. In particular, 
we define coeffsupp(O) :— 0, the empty set, and |coeffsupp(0)| :— —00. Also, we use 
coefF(P, k) Pk to denote the coefficient matrix (vector) Pk of z^ in P. In this paper, 
always denotes a general zero matrix whose size can be determined in the context. 

The Laurent polynomials that we shall consider in this paper have their coeffi- 
cients in a subfield F of the complex field C. Let F denote a subfield of C such that F 
is closed under the operations of complex conjugate of F and square roots of positive 
numbers in F. In other words, the subfield F of C satisfies the following properties: 

xeF and ^/ye¥, Vx,2/eF with y > 0. (1.6) 

Two particular examples of such subfields F are F = R (the field of real numbers) and 
F = C (the field of complex numbers). 

Now, we introduce the general matrix extension problem with symmetry. Through- 
out the paper, r and s denote two positive integers such that 1 ^ r ^ s. Let P be an 
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rxs matrix of Laurent polynomials with coefficients in F such that P(z)P* (z) ~ 1^ for 
all z e C\{0} and the symmetry of P is compatible, where denotes the r x r iden- 
tity matrix. The matrix extension problem with symmetry is to find an s x s square 
matrix P^ of Laurent polynomials with coefficients in F and with symmetry such that 
[/r.,0]Pe = P (that is, the submatrix of the first r rows of Pg is the given matrix P), 
the symmetry of Pg is compatible, and Pe(z)P*(z) = Is for all z e C\{0} (that is, Pg 
is paraunitary). Moreover, in many applications, it is often highly desirable that the 
coefficient support of Pe can be controlled by that of P in some way. 

In this paper, we study this general matrix extension problem with symmetry 
and we completely solve this problem as follows: 

Theorem 1. Let ¥ be a subfield of C such that (|1.6p holds. Let P be an r x s 
matrix of Laurent polynomials with coefficients in F such that the symmetry of P is 
compatible and P(z)P*(z) = for all z € C\{0}. Then there exists an s x s square 
matrix Pg , which can be constructed by Algorithm{l\in section 2 from the given matrix 
P, of Laurent polynomials with coefficients in F such that 

(i) [Ir, 0]Pe — P, that is, the submatrix of the first r rows of Pe is P; 

(ii) Pe is paraunitary: Pe(z)P*(z) — Is for all z G C\{0}; 

(iii) The symmetry of Pe is compatible; 

(iv) The coefficient support of Pe is controlled by that of P in the following sense: 

|coeffsupp([Pe]j,fc)| ^ max |coeffsupp([P]„,fc)|, l^j,k^s. (1.7) 

Theorem [T] on matrix extension with symmetry is built on a stronger result which 
represents any given paraunitary matrix having compatible symmetry by a simple 
cascade structure. The following result leads to a proof of Theorem [1] and completely 
characterizes any paraunitary matrix P in Theorem [1] 

Theorem 2. Let P be an r x s matrix of Laurent polynomials with coefficients in 
a subfield ¥ of C such that (|1.6p holds. Then P(z)P*(z) — Ir for all z G C\{0} and 
the symmetry of P is compatible as in ()1.4p , if and only if, there exist s x s matrices 
Pq, . . . , P.j+1 of Laurent polynomials with coefficients in ¥ such that 

(1) P can be represented as a product of Pq, . . . , Pj+i-' 

P(z) = [Ir,0]Pj+i{z)Pj{z) ■ ■ ■ Pi(z)Po(z); (1.8) 

(2) Pj, 1 ^ J ^ J are elementary: Pj{z)P*{z) — Ig and coeffsupp(Pj) C [—1, 1]; 

(3) (Pj+i, Pj) has mutually compatible symmetry for all ^ j ^ J; 

(4) Pq = UJg and Pj+i = diag{\}senls-r), where \^s9l, ^862 are products of a 
permutation matrix with a diagonal matrix of monomials, as defined in (j2.2p ; 

(5) max [|coeffsupp([P]m,rt)|/2] , where [•] is the ceiling function. 

The representation in (jl.Sp (without symmetry) is often called the cascade struc- 
ture in the literature of engineering, see [13l[T4j[T9]. In the context of wavelet analysis, 
matrix extension without symmetry has been discussed by Lawton, Lee and Shen in 
their interesting paper 15^ and a simple algorithm has been proposed there to de- 
rive a desired matrix Pe from a given row vector P of Laurent polynomials without 
symmetry. In electronic engineering, an algorithm using the cascade structure for 
matrix extension without symmetry has been given in |19| for filter banks with per- 
fect reconstruction property. The algorithms in jl5| 119) mainly deal with the special 
case that P is a row vector (that is, r = 1 in our case) without symmetry and the 
coefficient support of the derived matrix Pe indeed can be controlled by that of P. 
The algorithms in |15[ ITS] for the special case r = 1 can be employed to handle a 
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general r x s matrix P without symmetry, see [151 119) for detail. However, for the 
general case r > 1 , it is no longer clear whether the coefficient support of the derived 
matrix Pg obtained by the algorithms in [121 HH] can still be controlled by that of P. 

Several special cases of matrix extension with symmetry have been considered 
in the literature. For F = R and r = 1, matrix extension with symmetry has been 
considered in |16j . For r = 1, matrix extension with symmetry has been studied 
in [2 and a simple algorithm is given there. In the context of wavelet analysis, 
several particular cases of matrix extension with symmetry related to the construction 
of wavelets and multiwavelets have been investigated in [H [HI [Zl HI [SI IHl [HI [H] ■ 
However, for the general case of an r x s matrix, the approaches on matrix extension 
with symmetry in [71 [1^ for the particular case r = 1 cannot be employed to handle 
the general case. The algorithms in [71 [16] are very difficult to be generalized to the 
general case r > 1, partially due to the complicated relations of the symmetry patterns 
between different rows of P. For the general case of matrix extension with symmetry, 
it becomes much harder to control the coefficient support of the derived matrix Pg, 
comparing with the special case r — I. Extra effort is needed in any algorithm of 
deriving Pg so that its coefficient support can be controlled by that of P. 

The contributions of this paper lie in the following aspects. Firstly, we satisfac- 
torily solve the general matrix extension problem with symmetry for any r, s such 
that 1 ^ r ^ s. More importantly, we obtain a complete representation of any r x s 
paraunitary matrix P having compatible symmetry with 1 ^ r ^ s. This representa- 
tion leads to a step-by-step algorithm for deriving a desired matrix Pe from a given 
matrix P. Secondly, we obtain an optimal result in the sense of (|1.7p on controlling 
the coefficient support of the desired matrix Pg derived from a given matrix P by our 
algorithm. This is of importance in both theory and application, since short support 
of a filter or a multiwavelet is a highly desirable property and short support usually 
means a fast algorithm and simple implementation in practice. Thirdly, we introduce 
the notion of compatibility of symmetry, which plays a critical role in the study of 
the general matrix extension problem with symmetry for the multi-row case (r > 1). 
Fourthly, we provide a complete analysis and a systematic construction algorithm for 
d-band symmetric filter banks and symmetric orthonormal multiwavelets. Finally, 
most of the literature on the matrix extension problem only consider Laurent poly- 
nomials with coefficients in the special field C ([S]) or R {[IllIS]). In this paper, our 
setting is under a general field F, which can be any subfield of C satisfying (|1.6p . 

The structure of this paper is as follows. In section 2, we shall present a step-by- 
step algorithm which leads to constructive proofs of Theorems[T]and[2l In section 3, we 
shall discuss an application of our main results on matrix extension with symmetry to 
the design of symmetric filter banks in electronic engineering and to the construction of 
symmetric orthonormal multiwavelets in wavelet analysis. Examples will be provided 
to illustrate our algorithms. Finally, we shall prove Theorems [1] and [2] in section 4. 

2. An Algorithm for Matrix Extension with Symmetry. In this section, 
we present a step-by-step algorithm on matrix extension with symmetry to derive a 
desired matrix Pe in Theorem [2] from a given matrix P. Our algorithm has three steps: 
initialization, support reduction, and finalization. The step of initialization reduces 
the symmetry pattern of P to a standard form. The step of support reduction is the 
main body of the algorithm, producing a sequence of elementary matrices Ai, . . . , Aj 
that reduce the length of the coefficient support of P to 0. The step of finalization 
generates the desired matrix Pe as in Theorem [51 More precisely, our algorithm 
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written in the form of pseudo-code for Theorem [5] is as fohows: 

Algorithm 1. Input P as in Theorem,\^ with SP = {S9i)*Sd2 for some 1 x r 
and 1 X s row vectors 0i and 62 of Laurant polynomials with symmetry. 

1. Initiahzation: Let Q :— U^g^PU^gj- Then the symmetry pattern of Q is 

SQ — [Iri: Irgi^lrsj ^lr4] [isu ls27^ i ^ '^84]: (2-1) 

where all nonnegative integers ri, . . . , r4, si, . . . , S4 are uniquely determined by SP . 

2. Support Reduction: Let Pq '.— UJg,^ and J .= I. 
while (|coefFsupp(Q)| > 0) do %% outer while loop 

Let Qo Q, [fci,fc2] := coeffsupp(Q), and Aj Ig. 
if k2 — — fci then 

for j from 1 to r do 

Let q := [Qo]j.: and p :— [Q]j,:, i/ie jth rows of Qq and Q, respectively. 

Let [^1,^2] := coeffsupp(q), ^ := ^2 — ^1, and Bj Ig. 

if coeffsupp(q) = coeffsupp(p) and £ ^ 2 and {ii — ki or £2 = ^2) then 

Bj := Bq. Aj := AjBj. Qo := QqBj. 
end if 
end for 

Qo takes the form in (|2.8p . 

Let B(_fe2 fe^) := 7^, Qi := Qo, ji := 1 and j2 := ^3 + r4 + 1. 
while ji ^ ri + r2 and j2 ^ r do %% inner while loop 

Let qi := [Qi]j,^: and q2 ■= [Qi]j2,:- 

if coeff(qi, fci) — then ji ji + 1. end if 

if coeff(q2, fc2) — then j2 J2 + 1- end if 

if coeff(qi, fci) 7^ and coeff(q2, ^2) 7^ then 

B(-fe2,fe2) B(-fe2,fe2)B(qi,q2)- Ql := QlB(q,,q2). Aj := AjBjq^^q^). 

ji := + 1- i2 J2 + 1- 
end if 

end while XX end inner while loop 
end if 

Ql takes the form in (|2.8I) with either coeff(Qi, — fc) = or coeff(Qi, fc) = 0. 

Let Aj := AjBg^ and Q := QAj. 

Then SQ = [l^i , -lr2, zlrj, -^^Iril'^il^i , 

Replace si, . . . , S4 by s'^, . . . , s'^, respectively. Let P j :— Aj and J := J + 1. 
end while XX end outer while loop 

3. Finalization: Q = diag(i^i, 7^2, ^3, ^4) for some rj x Sj constant matrices Fj in ¥, 
j = 1,...,4. LetU := diag(C/Fi,?7j^2,f7i^3,f7F4) so that QU = [7^,0]. DefinePj t/* 
anrf P.7+1 dia.g{Usei, L-r)- 

Output a desired matrix Pg satisfying all the properties in Theorem\^ 

In the folfowing subsections, we present detailed constructions of the matrices 
U^e, Bq, B(qj q^), Bqj, aud Up appearing in Algorithm [TJ 

2.1. Initialization. Let be a 1 xn row vector of Laurent polynomials with sym- 
metry such that SO = [siz'^^ , . . . ,e„z^"] for some ei, . . . ,e„ € {—1, 1} and ci, . . . ,c„ S 
Z. Then, the symmetry of any entry in the vector Od\a,g(z~^^^^^^ , ■ ■ ■ , z~^'^"/^'^ ) be- 
longs to {±1,±2~^}. Thus, there is a permutation matrix Eg to regroup these four 
types of symmetries together so that 

S{e^se) - [I„i,-l„2,^"'ln3,-^"'ln4], (2.2) 
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where Ugg :— diag(z^rci/2l ^ _ _ ^ ^-\c„/2'] -j^^^ denotes the 1 x to row vector [1, . . . , 1], 
and ni, . . . , are nonnegative integers uniquely determined by S6. Since P satisfies 
(jl.4p . it is easy to see that Q := U^g^PU^e^ has the symmetry pattern as in (|2.ip . 
Note that Use^ and U^ga do not increase the length of the coefficient support of P. 



2.2. Support Reduction. Denote Q := UJg PU^^a in Algorithm [T] The 
outer while loop in the step of support reduction produces a sequence of elementary 
paraunitary matrices Ai , . . . , A j that reduce the length of the coefficient support of 
Q gradually to 0. The construction of each A^- has three parts: {Bi, . . . , 6^}, B(_/; 
and Bqj. The first part {Bi, . . . , B^} (see the for loop) is constructed recursively for 
each of the r rows of Q so that Qo '■— QBi • • • B^ has a special form as in (|2.8p . If 
both coeff(Qo,— fc) ^ and cocff(Qo,fc) ^ 0, then the second part B(^_i.,k) (see the 
inner while loop) is further constructed so that Qi := QqB^^j, takes the form in 
(|2.8p with at least one of coeff(Qi, — fc) and coeff(Qi,fc) being 0. Bq^ is constructed 
to handle the case that coeffsupp(Qi) — [— fc, fc — 1] or coeffsupp(Qi) = [—k + 1, k] so 
that coefl[supp(QiBQj C [-fc + l,fc- 1]. 

Let q denote an arbitrary row of Q with |coeffsupp(q)| ^ 2. We first explain how 
to construct Bq for a given row q such that Bq reduces the length of the coefficient 
support of q by 2 and keeps its symmetry pattern. Note that in the for loop, Bj is 
simply Bq with q being the current jth row of QBq ■ ■ ■ Bj_i, where Bq := Is- 

By (|2.ip . we have Sq = ez'^[lsi, ~ls2, z^^'^ss, —z^^'i-si] for some e G {—1,1} 
and c E {0, 1}. For e = —1, there is a permutation matrix such that S{qE^) = 
z'^[ls2, —isi, z~^ls4, —z^^'i-S'j]- For e = 1, we let := Is- Then, qE^ must take the 
form in either ([Q)) or ([2^ with f i ^ as follows: 



£2-2 

q£;e=[fi,-f2,gi,-g2]/^ + [f3,-f4,g3,-g4]/^+'+ coeS{qE,,e)z' 
+ [f 3, f 4, gi,g2W'''^ + [f 1, f 2, 0, 0]z^=; 

^2-2 

qS, =[0,0,fi,-f2]/^ + [gi,-g2,f3,-f4]/^+'+ J2 coeS{qE„£)z' 
+ [g3 , g4 , f 3 , f 4] " ' + [gi , g2 , f 1 , f 2] 2:'^ 



If qi?e takes the form in (12. 4p . we further construct a permutation matrix i?q such that 
[gi,g2,fi,f2]i^q = [fi,f2,gi,g2] and define Uq,e := E^Eqdiiig{Is-s^, z~^Is^), where 
Sg is the size of the row vector [gi, g2]- Then qUq^e takes the form in (|2.3p . For qE^ of 
form (|2.3p . we simply let Uq,e := E^. In this way, qo ■— qUq,e always takes the form 
in (1131) with f 1 7^ 0. ' 

Note that U^.e^* ^ = h and ||f i|| = ||f2|| if qoqS = 1, where ||f || := ^/fF. Now 
we construct an s x s paraunitary matrix Bq„ to reduce the coefficient support of qo 
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as m 



from [^1,^2] to [ii + 1,^2 - 1] as follows: 



B! := - 



cFi 





Rid + -) 










Cfj 

CF2 


-gl(l-i) 




-g2(l + i) 




|i^fi(l + z) 



-|^f2(l-^) 




cGi 






^fi(l ~ z) 



-|^f2(l + ^) 








Cg^g2 
CG2 



, (2.5) 



where Cf ^ 



Igil 



if gi + 0; 

c otherwise, 

2 , r,_2 , r,_2 , i_ 120/2 



g2||, Co := coeff(qo,^i + l)coeff(qo, -^2)/cfi, 

(2.6) 



2ef^ -CQ 
C 



if g2 / 0; 
otherwise, 



and 



1, 2 are unitary constant extension matri- 



^,^*] = C/f^-,[g;,G*] = ^g.forj 
ces in F for vectors f ^ , gj in F, respectively (see section 4 for a concrete construction 
of such unitary matrices \Ji- and C/g^ ). Here, the role of a unitary constant matrix \Jt 
in F is to reduce the number of nonzero entries in f such that fC/f — [||f||,0,...,0]. 
The operations for the emptyset are defined by|j0|j=0, + A = A and % ■ A — % 
for any object A. 

Define Bq := Uq.eBq„U* ^. Then Bq is paraunitary. Due to the particular form of 
Bqji as in (j2.5p . direct computations yield the following very important properties of 
the paraunitary matrix Bqi 



(PI) 5Bq = [1, 



"ls2 7 Z\g 



'^\sX coeffsupp(Bq) 



[—1,1], and coeffsupp(qBq) = [£i + 1, ^2 — 1]- That is, Bq has compatible sym- 
metry with coefficient support on [—1,1] and Bq reduces the length of the 
coefficient support of q exactly by 2. Moreover, 5(qBq) = 5q. 
(P2) if (p, q*) has mutually compatible symmetry and pq* = 0, then 5(pBq) = 5(p) 
and coeffsupp(pBq) C coeffsupp(p). That is, Bq keeps the symmetry pattern 
of p and does not increase the length of the coefficient support of p. 
Next, let us explain the construction of B(_j, ^.j. For coeffsupp(Q) = [— A:,fc] with 
fc ^ 1, Q is of the form as follows: 



-Fii 
--F12 



--F2: 
F22 



G31 
— G32 



G42 



F31 
-F32 



-F41 
F42 



--F52 



G71 
-G7 



— Gsi 

G82 



Gil 
— G12 



-G2 
G22 



F71 
-F72 



-Fgi 

F82 



+ ^ cocff(Q, 



+ 


F52 


Fei 
Fei 


G31 
G32 


G41 
G42 


fc-i , 

z + 


r Fii 
F12 


F21 
F22 










k 


G51 
. G52 


Gei 

G62 


F71 
F72 


Fsi 

Fs2 . 


Gil 
. G12 


G21 
G22 


F31 
F32 


F41 
F42 . 


z 



(2.7) 



with all FjkS and Gjk& being constant matrices in F and -Fn, -F22, -F42 being of 
size ri x si, r2 x S2, ^'3 x S3, ^4 x S4, respectively. Due to Property (PI) and (P2) of 
Bq, the for loop in Algorithm [1] reduces Q in (|2.7p to Qq := QBi • • • B^ as follows: 









G31 


— G41 








— G32 


G42 
































Gil 
G12 



G21 
G22 
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If either coeff(Qc 
and B 



I. 



-k) = or coefF(Qo, k) = 0, then the inner while loop does nothing 
If both coeff(Qo,— fc) 7^ and coeff(Qo,fc) 7^ 0, then Bi^-k,k) is 
constructed recursively from pairs (qi, q2) with qi, q2 being two rows of Qo satisfying 
coefF(qi,— fc) 7^ and coeff(q2,fc) 7^ 0. The construction of B(qj q^) with respect to 
such a pair (qi, q2) in the inner while loop is as follows. 

Similar to the discussion before (|2.3|) . there is a permutation matrix ^'(qi.qj) such 
that qi := qiE^^^ q^j and q2 := q2£'(qi.q2) take the following form: 



[ q2 J [000 

+ ^ cooff'( 



94 



-fe 



97 



esi -£92 e/7 
fe 93 94 



96 fr fs 



-98 



£/8 








31 92 



(2.9) 



where e € {—1, 1} and all gj's are nonzero row vectors. Note that ||gi|| = ||g2|| =: c 
and llgall = ||g4|| =: Cgg. Construct an s x s paraunitary matrix B(qj q^j as follows: 





SI 





g3(l + i) 


g4(l-i) " 




cGi 
















'^Sl 


-g3(l-i) 


-g4(l + i) 


1 













c 


3^gi(l + 2) 

El 


-^g2(l-^) 
SI 


-^g3 

S3 













cGa 







l^gl(l-^) 

SI 


-|^g2(l + ^) 
SI 





-^g4 
S3 













CG4 



(2.10) 



where cq := coeff(qi, -fc + l)coeff(q^, -fc)/cg, , c := (Icop + 4c|)i/2^ and [^,G*] = 
Ug- are unitary constant extension matrices in F for vectors gj in F, j = 1, ... ,4, 
respectively. Let B(q^^q2) := ii;(q^_q2)B(q^^q2)i?^^_q^). Similar to Property (PI) and 
(P2) of Bq, we have the following very important properties of B(qj q^): 

(P3) 5B(qj_q2) = [Isj,-ls2,zls3,-2lsj^[ls^,-ls2,z"^ls3,-z"^ls4], the Coeffi- 
cient support of B(qj_q2) is on [—1, 1], coeffsupp(qiB(qj^q2)) C [— fc + 1, fc — 1] 
and coeffsupp(q2B(qj qj.)) C [— fc + l,fc— 1]. That is, B(qj q^) has compatible 
symmetry with coefficient support on [—1, 1] and B(qj q,) reduces the length of 
both the coefficient supports of qi and q2 by 2. Moreover, 5(qiB(qj q^)) = Sqi 
and 5(q2B(q^^q,)) = 5q2. 
(P4) if both (p,q|) and (PjqJ) have mutually compatible symmetry and pq^; = 
pq2 = 0, then 5(pB(q^ q^)) = Sp and coeffsupp(pB(q^ q^j) C coeffsupp(p). 
That is, Bjq^qj) keeps the symmetry pattern of p and docs not increase the 
length of the coefficient support of p. 
Now, due to the Property (P3) and (P4) of B(qj q^), B(^_f; k) constructed in the 
inner while loop reduces Qo of the form in (|2.8[) with both coeff(Qo,— fc) 7^ and 
coeff(Qo,fc) 7^ 0, to Qi QQB(-k,k) of the form in (|2.8p with either coeff(Qi,— fc) = 
coeff(Qi,fc) — (for this case, simply let Bq^ := Is) or one of coeff(Qi,— fc) and 
coeff(Qi, fc) is nonzero. For the latter case, Bq^ := diag(C/iWi, /s3+s^)£^ with matrices 
[/i,Wi constructed with respect to coeff(Qi,fc) 7^ or Bq^ := diag(/sj+s2 , [/3W3)i? 
with C/3, W3 constructed with respect to coeff(Qi, — fc) 7^ 0, where i? is a permutation 
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matrix. Bq 
form in 



is constructed so that coeffsupp(QiBQ J C [—k + l,fc — 1]. Let Qi take 
. The matrices Ui,\Ni or J/a, W3, and E are constructed as follows. 



Let Ui :— diag(/7g^, U^J and U3 :— diag([/g^, UgJ with 



" Gil ■ 


, G2 


G21 


, G3 :— 


G31 


, G4 :— 


G41 


G12 




G22 




G32 




G42 



(2.11) 



Here, for a nonzero matrix G with rank m, Uq is a unitary matrix such that GUq = 
[R, 0] for some matrix R of rank m. For G = 0, C/g '■— I and for G = 0, [/g := 0. 
When GiG]; = G2G2, J7gi and fJcj can be constructed such that Gif/ci = 0] and 
G2UG2 = 0] (see section 4 for more detail). 

Let mi, TO3 be the ranks of Gi, G3, respectively (mi = when coefF(Qi,fc) = 
and ma = when coeff(Qi,-fc) = 0). Note that GiG^ = GaG^ or GgG^ = G4GI 
due to QiQ| = Ir- The matrices Wi, W3 are then constructed by: 





r Ui 




U2 






r U3 




U4 




Wi — 




is 1 — m 1 






,W3 := 










U2 




Ui 




U4 




Us 






















Isi — m^ _ 



(2.12) 



-U2(-z) := 



2 and 113(0) = U4(-2) := ^^Im^- 



where Ui(z) 

Let Wqj := diag([/iWi, /s3+s^) for the case that coeff(Qi,fc) ^ or Wq^ := 
diag(Jsi+s2, C/3W3) for the case that coeff(Qi, —k) ^ 0. Then Wq^ is paraunitary. By 
the symmetry pattern and orthogonality of Qi, Wq^ reduces the coefficient support of 
Qi to [— fc + 1, fc — 1], i.e., coeffsupp(QiWQ J = [— fc + 1, fc — 1]. Moreover, Wq^ changes 



mi 1 Ima i ■ 



the symmetry pattern of Qi such that iS(QiWqJ 
E is then the permutation matrix such that 

cS(QiWqJ£; = [1, 



[1. 



1-^2 7 ^ Irs ; ^lr4 



with so 



•-Si —mi +7713 ' • 



■-52 — 77li+m3 J 



■-Ss — 7773+7711 



1^2 1 , ■ 

-7713 ; lm3 1 " 

1-1 

-■-S4 — 7773 + 7771 



-zlr^^SBi with 



{sei)E. 



3. Application to Filter Banks and Orthonormal Multiv^ravelets with 
Symmetry. In this section, we shall discuss the application of our results on matrix 
extension with symmetry to d-band symmetric paraunitary filter banks in electronic 
engineering and to orthonormal multiwavelets with symmetry in wavelet analysis. In 
order to do so, let us introduce some definitions first. 

We say that d is a dilation factor if d is an integer with |d| > 1. Throughout this 
section, d denotes a dilation factor. For simplicity of presentation, we further assume 
that d is positive, while multiwavelets and filter banks with a negative dilation factor 
can be handled similarly by a slight modification of the statements in this paper. 

Let F be a subfield of C such that (0 holds. A low-pass fiher oq : Z ^^ f^'' 
with multiplicity r is a finitely supported sequence of r x r matrices on Z. The symbol 
of the filter oq is defined to be siq{z) := 'Ylik£z'^o{k)z^ ^ which is a matrix of Laurent 
polynomials with coefficients in F. Moreover, the d-band subsymbols of oq are defined 
by 80:7(2;) := Vd X^fcez ^^0(7 + dk)z'^, 7 G Z. We say that ao (or ao) is a d-band 
orthogonal filter if 



d-l 

^ao-^{z)a*o.^{z) = ^ 

7=0 
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z e C\{0}. 



(3.1) 



To design an orthogonal filter bank with the perfect reconstruction property, one has 
to design high-pass filters ai, . . . , a^^i : Z fr'^r g^g}j ^j-^g^^ ^]-jg polyphase matrix 



ao;o(-z) ••■ ao;d-i(2;) 
ai;o(-z) ••■ ai;d-i(z) 

ad-i;o(2) ••• ad-i;d-i(z) 



(3.2) 



is paraunitary, that is, V{z)V*{z) = I^r, where each am-.-y is a subsymbol of am for 
TO, 7 = 0, . . . , d — 1, respectively. Symmetry of the filters in a filter bank is a very 
much desirable property in many applications. We say that the low-pass filter ag (or 
flo) has symmetry if 

ao(z) = diag(eiz^=i , . . . , e,z^^'-)ao(l/z)diag(eiz-^S . . . , ErZ'^^) (3.3) 

for some £i, . . . , G {^Ij 1} and ci, . . . , G M such that dci — Cj G Z for all 
£,j = l,...,r. To design a symmetric filter bank with the perfect reconstruction 
property, from a given d-band orthogonal low-pass filter ao, one has to construct 
high-pass filters ai, . . . , a^-i : Z i— >■ F''^'' such that all of them have symmetry that is 
compatible with the symmetry of ao in (|3.3p and the polyphase matrix V in (|3.2p is 
paraunitary. 

For / G -Li(M), the Fourier transform used in this paper is defined to be /(^) := 
J-g^f{x)e~^^^dx and can be naturally extended to £2(K) functions. For a d-band 
orthogonal low-pass filter ao, we assume that there exists an orthogonal d-refinable 
function vector (j) = [ipi, . . . , ipr]'^ associated with the low-pass filter ao, with compactly 
supported functions 0i, . . . , 1/)^ in i2(]R), such that 

0(dO = ao(e"^«)^(C), with ||0(O)|| = 1, (3.4) 

and 



((/)(•- fc),0) := (j){x - k)4>ix) dx = S{k)Ir, keZ, (3.5) 

where 6 denotes the Dirac sequence such that (5(0) = 1 and 5{k) = for all k ^ Q. 
Define multiwavelet function vectors -0™ = VPT i • ■ • ; i'TV associated with the high- 
pass filters a„i, TO = 1, . . . , d — 1, by 

V^(dC) := a„(e-'«)0(e), ^ m = 1, . . . ,d ~ 1. (3.6) 

It is well known that {V'^i ■ • ■ generates an orthonormal multiwavelet basis in 

L2(M); that is, {d^l^i^'^id^ ■ -k) : j. A: G Z; m = 1, . . . , d - 1; £ = 1, . . . , r} is an 
orthonormal basis of L2(M), for example, see [H [3 [11] [17] and references therein. 

If ao has symmetry as in p.3p and if 1 is a simple eigenvalue of ao(l), then it is 
well known that the d-refinable function vector (j) in (|3.4p associated with the low-pass 
filter ao has the following symmetry: 

0l(ci - •) = ^l*/*!! 02(C2 - •) = ^^202, (l)r{Cr - ■) = er(t)r- (3.7) 

Under the symmetry condition in (j3.3|) . to apply Theorem [l] we first show that 
there exists a suitable paraunitary matrix U acting on Pa„ := [ao;Oi ■ ■ ■ , ao;d-i] so that 
Pan U has compatible symmetry. Note that Rao itself may not have any symmetry. 

Lemma 1. Let Pag := [ao;o, ■ • ■ j ao;d-i], where ao;o, • ■ • , ao;d-i are d-band sub- 
symbols of a d-band orthogonal filter ao satisfying p.3p . Then there exists a dr x dr 
paraunitary matrix U such that PaoU has compatible symmetry. 
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Proof. From p.3p . we deduce that 
[ao;7(z)]f,i = eiejZ^l:>[ao.QjXz~^)]ej, 7 = 0, . . . , d - 1; f , j = 1, . . . , r, (3.8) 
where 7, Qj^- G F {0, . . . , d — 1} and R'J j, Qj j are uniquely determined by 

dee - - 7 = dRl^ + Ql^ with £ Z, Qj^. G F. (3.9) 

Since dce — Cj G Z for all j = 1, . . . , r, we have ce — cj G Z for all ^, = 1, . . . , r 



and therefore, ^ is independent of ^. Consequently, by p.Sp . for every 1 ^ j ^ r, 
the jth column of the matrix ao;-y is a flipped version of the jth column of the matrix 
ao;Qj ■■ Let Kj^j G Z be an integer such that |coeffsupp([ao;'y]: j + z^J-"' [ao;Qj is 
as small as possible. Define P := [bo;o, ■ • ■ , bo;d-i] as follows: 

{[aO;7]:,i: 7 = Q]/, 

73([aO;7]:j+^"^'iao;Q7j:,,), 1< Qlp (3.10) 
71 ([aO;7]:j - ^''^■■iao;Qj_J:,j), 7 > Q],j, 

where [ao;^]:^ denotes the jth column of ao;^. Let U denote the unique transform 
matrix corresponding to (I3.10p such that P := [bo;o, • • . , bo;d-i] = [ao;Oj • • • 7 ao;ci-i]U. 
It is evident that U is paraunitary and P = PaoU. We now show that P has compatible 
symmetry. Indeed, by (I3.8p and (|3.10l) . 

[Sho,^]L, = sgn{Ql^ - 7)£,e,z«L+«-^, (3.11) 

where sgn(a;) = 1 for x ^ and sgn(2;) = —1 for a; < 0. By p.9p and noting that Qjj 
is independent of £, we have 

for all 1 ^ n ^ r, which is equivalent to saying that P has compatible symmetry. □ 
Now, for a d-band orthogonal low-pass filter ao satisfying p.3p . we have the 
following algorithm to construct high-pass filters ai, . . . , ad-i such that they form a 
symmetric paraunitary filter bank with the perfect reconstruction property. 

Algorithm 2. Input an orthogonal d-band filter Sq with symmetry in (j3.3p . 

(1) Construct U as in p.lOp such that P := PaoU 'f^-s compatible symmetry: 
SP — [eiz'^i , . . . , e^z'^'^J^iS^ for some ki,...,kr G Z and some 1 x dr row 
vector 9 of Laurent polynomials with symmetry. 

(2) Derive Pg with all the properties as in Theorem{l\ from P by Algorithm{l\ 

(3) LetV :— PgU* —: (a,„;^)o^m,7^d-i as p.2p . Define high-pass filters 

1 

a™(z) := — Va„^^(z^)z^, m = l,...,d-l. (3.12) 

Output a symmetric filter bank {ao, ai, . . . , ad-i} with the perfect reconstruction prop- 
erty, i.e. V in (j3.2l) is paraunitary and all filters a™, m = 1, . . . ,d — 1, have symmetry: 

a™(z) = diag(£rz^='", . . . ,e;!'z^<")a„(l/z) diag(eiz-^\ . . .,erZ'^^), (3.13) 
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where c^' 
rn = 1, . 



(fc™ - fc^) + Q e K and all e {-1, 1}, fc]," e Z, /or l = \,. 
~ 1, are determined by the symmetry pattern of Pg as follows: 



[eiz 



ki 



, . . . , Or 



1 k'^ 



]^56l := SPe 



, r and 



(3.14) 



Proof. Rewrite Pe — (bm;-y)o<5m,7<d-i as a d x d block matrix with r x r 
blocks bm;7- Since Pe has compatible symmetry as in (|3.14l) . we have [Sbm;'y]e,: = 
£fetz^'^-'^'[SbQ,^]i. for £ = 1, . . . ,r and m = 1, . . . ,d - 1. By (|3TT|l . we have 



By (pTTSl and the definition of U* in (|3?T0l) . we deduce that 



(3.15) 
(3.16) 



This implies that [5a,„]£j- = £™£^.;2d(fej"-fef+cf)-cj^ which is equivalent to (|3A3)) with 
c™ := fc™ — fcf + C£ for TO = 1, . . . , d — 1 and £ = 1, . . . ,r. □ 

Since the high-pass filters ai, . . . , ad-i satisfy p.lSp . it is easy to verify that each 
ip"^ = . . . , V'™]^ defined in p.6p also has the following symmetry: 



^5"(cf -•) ^erv^r. V'r(c"--) = erv': 



2 J 



In the following, let us present several examples to demonstrate our results and 
illustrate our algorithms. 

Example 1. Let d = 2 and r = 2. A 2-band orthogonal low-pass filter ao with 
multiplicity 2 in [S] is given by 



1 

40 



12(1 + ^-1) 16\/2z-i 
-V^(z2-9z-9 + z-i) -2(3z-10 + 3z-i) 



The low-pass filter ao satisfies p.3p with ci — — 1,C2 ~ and ei — 82 — 1- Using 
LemmalU we obtain Pa„ [ao;Oiao;i] and U as follows: 



P — _ 

^ 20 



6^/2 
9 - z IOV2 



6s/2 



-3\/2(l - 



' 1 





1 
















z 





—z 





_ 








^z 



Then P := PaoU satisfies 5P ~ [1, 2:]^[1, z ^, —1, 1] and is given by 

p _ \/2 r 6^/2 8\/2 
~20[4(l + z) 10 5(1 -z) -3(1 + z) 

Applying Algorithm [U we obtain a desired paraunitary matrix Pg as follows: 

6%/2 8\/2 

\/2 4(1 + z) 10 5(1 -z) -3(1 + z) 

' 20 4(1 + z) -10 5(1 -z) -3(1 + z) 

_ 4V^(1 - z) 5V^(z + 1) 3V2(z - 1) 

We have 5Pe = [1, z, z, — z]-^[l, z~^, — 1, 1] and coeffsupp([Pe]: j) C coeffsupp([P]. j) 
for all 1 ^ j ^ 4. Now, from the polyphase matrix V := PgU* =: (am;-y)o^m^^^i, we 
derive a high-pass filter ai as follows: 



ai(z) 



1 

40 



-V2(z2-9z-9 + z-i) -2(3z + 10 + 3z-i) 



2(z2-9z + 9-z-i) 
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Then the high-pass fiher ai satisfies (|3.13l) with cj 



and el — l,e\ = — 1. 



Example 2. Let d = 3 and r = 2. A 3-band orthogonal low-pass filter ao with 
multiplicity 2 in is given by: 



an (z) -I- flu (z 1) ai2{z) + z ^ai2{z 
a2i(z) -I- z'^a2i{z^'^) a22{z) + z'^a22{z^^) 



where 



an (2) = 90 + (55 - 5741)2 - (8 + 2741)2^ + (7^ - 47)2*; 
012(2) = 145 + 5\/4T + (1 - %/4T)2^ + (34 - 4^41)2^ 
a2i(2) = (111 + 9V4T)2^ + (69 - 9\/4T)2*; 
a22(2) = 9O2 + (63 - 3\/4T)2^ + (3\/4T - 63)2^ 

The low-pass filter ag satisfies (I3.3P with ci = 0,C2 = 1 and ei = £2 
Pao [ao;o, ao;i, ao;2], the matrix U constructed by Lemma[T]is given by 

72 





1 . From 



U := 



1 

71 





1 
1 
V2 
z 
z 





1 


— z 






1 





— z 



Let 

Co = 11 
ii6 = 3co; 
t53 = 400\/6/co; 



ti2 = 5(7- V41); 
^l5 = 3(3^/4T-13); 



C12 = 10(29 + V41); 
^25 = 6(7 + 3\/4T); 
f55 = 12\/6(\/4T- 1); t56 = 6\/6(4 + \/4T); 



tl3 
^26 
C66 



-5co; 

6(21-^/4T); 
3\/6(3 + 7V4T). 



Then P :— U satisfies SP — [1, z]-^[l, 1, 1, z ^, — 1, — 1] and is given by 



V6 
1080 



180^2 6i2(z) 6i3(z) ii5(z-z-i) ti6(z-z-i) 

180(H-z) I8OV2 t25(l-z) t2Q{l-z) 



where 6i2(z) = ti2(z-f z ^) + ci2 and hiz{z) — ti3(z — 2- 
we obtain a desired paraunitary matrix Pg as follows: 



Applying AlgorithmfT] 





" I8OV2 



612(2) 



613(2) 

180(H-2) 




180%/2 


il5(2- 1) 
t25(l-2) 


tie (2 — 

t26(l- 


1) 1 

^) 


1080 


360 



i'12(z) 




fcl3(z) 
V2 

90^(1 + 2) 




-360 




72 U 
*26 (-1 _ 
72 


^) 
^) 








V&tvi(l - 2) 

^(i-2) 


t53(l-2) 
^(i-2) 






t55(l + 2) 
665(2) 


t56(l + 
666(2 


^) 



where b^^iz) ~ — \/6(5ti5(z + z^^) + 3ci2)/10 and fe66(^) = —VGti(i{z + z^^)/2 + cge- 
Note that iSPe = [1, z, 1, z, — z, — 1]^[1, 1, 1, z~^, —1, —1] and the coefficient support of 
Pe satisfies coefI'supp([Pe]:.j) C cocffsupp([P]:^j) for all 1 ^ j ^ 6. From the polyphase 
matrix V := PgU* =: (am;-y)o^,„^^^2, we derive two high-pass filters ai,a2 as follows: 

V2 \ an(^) +a}i(z-i) al^iz) + z-^al^{z'^) 

1080 [ al,{z) + z^al.iz-^) al,{z) + z^a\2{z-') 

%/6 \ o\^{z)- z^al^{z-^) a\2{z)~ z^a\2{z-^) 

1080 [ al,{z) - ali(z-i) a^z) z'^a^z'^) 
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ai(z) 
^i{z) 



where 



ail (z) = (47 - 7\/4T)/ + 2(4 + \/4T)z^ + 5(\/4T - ll)z + 180; 
al^iz) = 2(2\/4T - 17)z^ + {V^l - 1)2^ - 5(29 + 

aliiz) = 3(37 + 3\/4T)z + 3(23 - 3Vil)z~^; 

^22(2) = -ISOz + 3(21 - V41) - 3(21 - 741)2"^ 

an(z) = (43 + 17^41)2 + (67 - 7^)^"^ 

alziz) = 11\/4T-31- (79 + \/4T)z"\ 

a^i(z) = (47 - 7^41)2* + 2(4 + V^)z'^ - 3(29 + V^)z- 

aliiz) = 2(2%/4T - 17)2^ + {V41 - 1)2^ + 3(3 + 7\/4T). 

Then the high-pass filters ai, 32 satisfy p.l3p with c\ = 0,C2 — I, £1 — £2 — ^ ^'^d 
= 1, c| = 0, = e| = — 1. 

As demonstrated by the following example, our Algorithm [2] also applies to low- 
pass filters with symmetry patterns other than those in p.3p . 

Example 3. Let d = 3 and r = 2. A 3-band orthogonal low-pass filter ap with 
multiplicity 2 in |5] is given by 

011(2;) 012(2:) 
021(2) 022(2:) J ' 

where 

011(2) = (11 - 14^/T7)2^ + (29 + 8\/T7)2 + 234 + (85 - 16\/T7)2"^ - (17 + 2\/T7)2"^ 
012(2) = (bVri - 16)2^ + {2 + \/T7)2^ + 238 - llVTl + (136 + 29\/T7)2"^ 
021(2) = (136 + 29\/T7)2^ + (238 - llVr7)z + {2 + VT7)z''^ + (dVTl - 16)2"^ 
022(2) = (-17 - 2\/T7)2^ + (85 - 16\/T7)2^ + 2342 + 29 + sVTl + (11 - 14\/T7)2"\ 

This low- pass filter ao does not satisfy (I3.3p . However, we can employ a very simple 
orthogonal transform := ^[ I _\ ] to ao so that the symmetry in (|3.3p holds. 
That is, for 30(2;) '■— EaQ{z)E, it is easy to verify that Sq satisfies (13.31) with ci = 
C2 = 1/2 and £1 = l,£2 = —1. Construct P^q := [a^0;07 a^o;i, a^0;2] from ao. The matrix 
U constructed by Lemma [T] from P^^ is given by: 

-1 " 

2-1 

1 0-1000 
1 0-100- 
^/20 

V2 _ 

Then P := PsqU satisfies SP = — 2-i]^[l, —1, —1, 1, 1, —1] and is given by 




234(1 4-^) ti2(l-i) ^13(1-7) 117^/2(l + i) ti6(l - i) 

t2l(l-|) t22(l + |) t23(l + i) t24(l-i) t2s{l - ^) t26(l + i) 
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''^'^ ^ 702 



V2 



where c = and tjk^s are constants defined as foUows: 

ti2 = 3(11 - 7l7); ti3 = 3(717-89); fi6 = 15^2(2 + 7l7); 

tai = 13(717- 17); t22 = 6(2 + 717); t23 = 6(37 - 717); 

t24 = -13(1 + 717); t25 = -1372(8 + 717); t26 = -372(7 + 107T7). 

Applying Algorithm [T] to P, we obtain a desired paraunitary matrix Pg as follows: 





" 234(1 + i) 
fei(l-i) 


tl2(l-i) 
t22(l + i) 


t23(l+i) 




i24(l-i) 


11772(1 + i) 

t25(l-i) 


ii6(l-7) ' 

t26(l+i) 


Pe = C 


^3i(l-|) 

i4l(l+i) 


t32(l + i) 
i42(l-i) 


t33(l+|) 
i43(l-^) 


i34(l-|) 
t44(l + i) 


t35(l-|) 
-72t4l(l+i) 


t36(l+|) 
i46(l-i) 









t62 




^63 


— 273t4i 



"^*44 






^66 



where all tjk 's are constants given by: 







726(61 + 25717) /4; 


t32 


= -3726(397 + 237T7)/52; 


^33 


= 3^ 


/26(553 + 237T7)/52; 


^34 


= 25726(1 + 717) /4; 


t35 


= 7 


13(25717 - 43)/2; 


^36 


= 15713(23717- 19)/26 


til 


= 9^ 


/26(1 - 3717) /4; 


^42 


= -3726(383 + 297T7)/52; 


t43 


= 3^ 


/26(297T7 + 227)/52; 


t44 


= 27726(1 + 717) /4; 


t46 


= 3^ 


/l3(1457T7-61)/26; 


^62 


= 9778(4l7l7-9)/26; 


^63 


= 9^ 


/78(ll7T7 + 9)/26; 


t(>6 


= 2773(717 + 15)/7T3. 



Note that Pe satisfies SP^ = z^^, 1, -1, -1, 1, 1, -1] and we 

have coefFsupp([Pe]: j) C coeffsupp([P]:.j) for all 1 ^ j ^ 6. From the polyphase 
matrix V :— PeU*, we derive two high-pass filters ai,af2 as follows: 



ai(z) 
32 (z) 



26 



36504 
778 



4056 



aliiz) + za^i(z-i) al^iz) - za^^iz^'^) 
' (z) al^iz) 



a 



ali(z) al^iz) 



where 



011(2) 


= (433 - 128717)2^ 


+ 13(25V 


17- 


43)2^ - (1226 + 197717)2 


012(2) 


= (128717-433)2^ 


+ 15(23V 


17- 


19)2^ - (758 + 197717)2; 


021(2) 


= 3(133 - 44717)2^ 


+ 117(3V 


Tr- 


1)2^ -3(73717 + 94)2; 


022(2) 


= 3(44717-133)2^ 


+ 3(145V 


17- 


61)2^ - 3(250 + 73717)2; 


0?i(2) 


= 13(1 + 717) (2=*- 


22^+2); 






012(2) 


= 13(3717- 1)(2^ - 


-2); 






2 1 \ 
021(2) 


= (9+ ll7T7)(2^ - 


2); 






2 / \ 
022(2) 


= (4l7T7-9)(2^ + 


242V137 


+ 18 


\/T72Vl37 + 2). 
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Then the high-pass fihers a^i and 32 satisfy p.l3p with cj = C2 = 1/2, e\ ~ —1, £2 — 1 
and c\= c^= 3/2, e} = l,e\ = —1, respectively. 

Let 31,32 be two high-pass fihers constructed from ^i,1f2 by 3i(z) Eai{z)E 
and 32(2) := Ea2{z)E. Then due to the orthogonahty of E, {30,31,32} stih forms 
a d-band fiher bank with the perfect reconstruction property but their symmetry 
patterns are different to those of ifo , a 1 1 sf2 ■ 

4. Proofs of Theorems [T] and [2l In this section, we shah prove Theorems 
[T] and [2] The key ingredient is to prove that the coefficient supports of Ai , . . . , A j 
constructed in Algorithm [T] are all contained inside [—1, 1]. Note that each Aj takes 
the form Aj = (Bi • • • Br)B(_fe fc) Bq^ . We first show that the coefficient support of 
B (Bi • • • B.r)B(_j. J.) is contained inside [—1, 1] and then show that the coefficient 
support of BBqj is also contained inside [—1, 1]. 

Let us first present a detailed construction for the unitary matrices Uf and Ug 
that are used in Algorithm [T] For a 1 x n row vector f in F such that ||f || 7^ 0, we 
define rif to be the number of nonzero entries in f and e^ :— [0, . . . , 0, 1, 0, . . . , 0] to 
be the jth unit coordinate row vector in R". Let Et be a permutation matrix such 
that fEf = [/i, . . . , 0, . . . , 0] with fj 7^ for j = 1, ... , rif . We define 

f In, if n-f 1; 

J^fr (4-1) 
I 1/1 1 



In - TK-WjVfVi , if > 1, 



where Vf :— f — |^||f ||ei. Observing that \\vf\\^ ~ 2||f||(||f|| — |/i|), we can verify 
that VtVf ~ In and fEtVi — ||f ||ei. Let Ui EfVt. Then C/f is unitary and satisfies 
Ut = for some (n - 1) x n matrix in F such that fC/f = [|lf ||, 0, . . . , 0]. 

We also define Uf := /„ if f = and C/f := if f = 0. Here, Uf plays the role of 
reducing the number of nonzero entries in f . More generally, for an r x n nonzero 
matrix G of rank m in F, employing the above procedure to each row of G, we can 
obtain an n x n unitary matrix Ug such that GUg — [R, 0] for some r x m lower 
triangular matrix R of rank m. If GiG^ = G2G2, then the above procedure produces 
two matrices Ug^Ug^ such that GiUgi — 0] and G2UG2 — [R,0] for some lower 
triangular matrix R of full rank. It is important to notice that the constructions of Uf 
and Ug only involve the nonzero entries of f and nonzero columns of G, respectively. 
In other words, we have 



[C/f],.: - = e,, if [f], = 0, 

pGh, ^ {[UgI.V ^ e„ if[G],, -0. 



(4.2) 



Next, we establish the following lemma, which is needed later to show that the 
coefficient support of (Bi • • • Br)B(_j,jj) is contained inside [—1, 1]. 

Lemma 1. Suppose B is an s x s paraunitary matrix such that coeffsupp(B) C 
[—1, 1] and SB = {S9)*S9 with SO = [Is^ , — Isa , "'^Iss i ^^^^154] for some nonnega- 
tive integers si, . . . , S4 such that si + S2 + S3 + s^ = s. Then the following statements 
hold. 

(1) Let p be a 1 X s row vector of Laurent polynomials with symmetry such that 
pp* — I, coeffsupp(p) = [fci,A;2] with k2 ~ ki ^ 2, and Sp — ez'^SO for some 
£ e {—1,1} and c e {0,1}. Let q := pB. // coeffsupp(q) — coeffsupp(p), 
then coeffsupp(BBq) C [—1,1], where Bq is constructed with respect to q as 
in section 2. 

(2) Let pi,P2 be two 1 x s row vectors of Laurent polynomials with symmetry 
such that PjiP*2 = 5{ji -32) for ji,j2 = 1,2, Spi = eiS9 and Sp2 = £2^56* 
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for some £1,62 G { — 1,1}, and coeffsupp(pi) = coeffsupp(p2) C [— fc, fc] with 
k ^ 1. Let qi piB and q2 '■— P2B. // coefFsupp(qi) — [— fc, fc— 1] and 
coeffsupp(q2) = [—k + l,k], i/ien coeffsupp(BB(qj q^)) C [—1,1], where B(q^ 
is constructed with respect to the pair (qi,q2) as in section 2. 
Proof. Due to Sp — ez^SO, as we discussed in section 2, there is an Up^^ such that 
pUp^e takes the form in (|2.3p . Since Up^e is a product of a permutation matrix and 
a diagonal matrix of monomials, we shall consider the case that Up^e = Is, while the 
proofs for other cases of Up.e can be obtained accordingly. Then p takes the standard 
form in ([Q]) with fi ^ 0. In this case, si > and S2 > due to ||fi|| = |lf2|| ^ 0. 
By our assumptions, q := pB must take the following form: 

fc2-2 

q := pB =[fi, -?2,gi, -g2]2'°' + [fa, -f4,g3, -g4]2''i'''^ + ^ coeff(pB, n)2:" 

n—k-i+2 

'=2-1 I n nU*2 



+ [f3,f4,gl,g2]2: 



with fl ^ 0. Then Bq is given by dp]) withjfi, ±2, gi, g2, i^i, F2, Gi, G2 
being replaced by fi, f2, gi, g2, ^i, F2, Gi, G2 respectively and all constants 
"^f 1 ' "^gi ' "^^2 : Co, c, Cg'^ , Cg^ being defined accordingly. 

Also, due to the symmetry pattern and coeffsupp(B) C [—1, 1], B is of the form: 







B3(l + i) 


B4(l-i) 1 




Ai{z+^) + D2 


C3(l-i) 


C4(l+i) 


B^{l + z) 


Cl(l-2) 




MZ-^,) 


£2(1-2) 


C2(l + 2) 




As{z+\) + Di _ 



(4.3) 



X Si 



where Aj's, -Bj's, Cj's and -Dj's are all constant matrices in F and Dj is of size Sj 
for j = 1,...,4. 

Let X ;= {1, si + 1, (1 - (5(s3))(si + S2 + 1), (1 - 5{si)){si + S2 + S3 + 1)} be 
an index set. It is easy to verify that coeffsupp([BBq]: j) C [—1,1] for all j ^ I. 
Hence, by coeffsupp(BBq) C [—2,2], we only need to compute coeff([BBq]: j, 2) and 
coeff([BBq]: J, — 2) for those j G I. Let us show that coeff([BBq]: j-, 2) = for j = 1, 
i.e., the coefficient vector of for the first column of BBq is 0. By coeff(pB,fci) = 
coeff(p, fci + l)coeff(B, —1) + coeff(p, fci)coeff(B, 0), we have 



f 1 = f 3A1 + tiA2 + tiDi + giBi - g2B2; 

f 2 = f 3^3 + tiAi + f 2D2 - glCl + g2C2; 

gl = f'iB'i + fiCs + gsAs + g4A6 + f 1B3 - f 2C3 + giD3; 

g2 = f 3-B4 + f 4C4 + g3^7 + g4As - f I-B4 + f 2C4 + g2-D4. 

Similarly, by coeff(BBq, 2) = coefF(B, l)coefF(Bq, 1), we have 



(4.4) 



off([BB, 



.2) 



A2 
Bi 
-S2 



A3 
Ai 
-Ci 
C2 






As 
Ae 






A^ 
An 











-f * 


1 




. -82 . 


c 



Alf * - A3I; 
A2fl - A4I* 
BifJ + Cif| + Asgj; - A^^ 
-Baf J - Cif ^ + Aeel - Assi 



Due to BB* = Is, we obtain 



AiAl - ^3^3 = 0, A1A2 - A3AI = 0; 
AiDt + DiAl + B'iBt - B4BI = 0; 
DiA*2 - AaD^ + B3C3* - B4C4* = 0; 
AiBt + AaCt + B3AI - BiA*j = 0; 
-AiB^2 - AaCi + B3AI - B4AI = 0. 
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Applying the above identities to ^if ^ — 2 and using (|4.4p . we get 

AJl - ^3?; = Ai{±3Ai + fiAz + ±iDi + giBi - g2-B2)* 

- ^3(13^3 + f 4^14 + f 2O2 - glCi + g2C2)"' 

= {AiAl - AsADn + {AiAl - AsAl)fl + {AiDl)fl 
+ {-A3D2)f; + (AiBl + A3C;)gl - {AiB; + A:iC;)g*2 

= -{DiAl + B3BI - BABl)fl - {DiAl + B3C3 - B4C4*)f2 

- (B3A5 - B4A*r)gl - {BsA; - BiA'8)g2 

= -Dl(f lAi + fiA^y - B3(f 1B3 + f2C3 + glAs + g2^6)* 
+ B4(f 1B4 + f 2C4 + gl^7 + g2^8)* = 0, 

where the last above identity follows by coeff(pB,fc2 + 1) = coeff(pB,A;i — 1) = 0. 
Similarly, we can show that ^2? | - ^4? 2 = 0: Bif^ + Cif 2 + ^5g^ - Ayg^ = 0, and 
-B2fl - Cif*2 + Aogl - Aszl = 0- Hence, coeff([BBq]:,i, 2) = 0. By similar computa- 
tions as above and using the paraunitary property of B, we have coeff([BBq]: j, ±2) = 
for all j € X. Therefore, we conclude that coeffsupp(BBq) C [—1, 1]. Item (1) holds. 

For item (2), up to a permutation matrix -E(qi,q2) as in section 2, B(qj^_q2) takes the 
form in p.lOp . Since B takes the form in (j4.3p . to show that the coefficient support 
of BB(_j, J,) is contained inside [—1, 1], we need to show that all the coefficient vectors 
^ig^ - ^3g2, Mgl - ^4g2, ^5g3 - ^7gl, ?aid A^gl - AggJ are zero. Again, using the 
paraunitary property of B and expressing gi , g2 , ga , g4 in terms of the original vectors 
from pi, p2 similar to (14. 4p . we conclude that coeffsupp(BB(qj q^)) C [—1, 1]. □ 

With the result of Lemma[I] the next lemma shows that the coefficient support of 
B := (Bi • • • Br)B(_fc is contained inside [—1, 1]. Moreover, the next lemma shows 
that the coefficient support of A := BBq^ is also contained inside [—1, 1]. 

Lemma 2. Suppose Q is anr x s matrix of Laurent polynomials such that QQ* = 
Ir, SQ satisfies (|2.ip . and coeffsupp(Q) — [^1,^2] with k2 — ki ^ 1. Then there exists 
an s X s paraunitary matrix A of Laurent polynomials with symmetry such that 

(1) coeffsupp(A) C [—1,1] and |coeffsupp(QA) I < |coeffsupp(Q)| — |coeffsupp(A)|; 

(2) if the jth column p :— [Q]:j of Q satisfies coeff(p, ki) = coeff(p, k2) = 0, then 
[A]j.: = = Gj . That is, any entry in the jth row or jth column of A 
is zero except that the (j, j) -entry [A]jj = 1; 

(3) SA = [ls^,-ls2,zlss,-zls^]'^[ls'^,-ls'^,z~^ls'^,-z~^ls'J for some nonneg- 
ative integers s'^, . . . , S4 such that s'l + + S3 + S4 = s. 

Proof. Let A — (Bi • • • Br)B(_j, j.) Bq^ be constructed as in Algorithm [TJ where 
Qi := Q(Bi • • • Br)B(_j, fc), B(_j, j.) is constructed in the inner while loop of Algorithm 
[1] and Bi, . . . , Bj. is constructed in the for loop of Algorithm [T] If /c2 7^ — fci, then 
Bi = • • • = Br = B(_fc fc) — Is and A is simply Bq^ , where Qi = Q is of the form in p.Sp 
with either coefF(Qi, — fc) = or coefF(Qi, k) ~ 0. In this case, by the construction 
of Bqj as in section 2, all items in Lemma [2] hold. We are already done. So, without 
loss of generality, we assume that k2 = —ki = k. 

We first show that the coefficient support of Bi • • • B^ is contained inside [—1, 1]. 
Let pj :— [Q]j.:, Bo :— Is, and :— PjBq ■ ■ ■ Bj_i for j = 1, . . . , r. Suppose we already 
show that coeffsupp(Bo • • • Bj_i) C [—1, 1] for j ^ 1. Then, according to Algorithm[l] 
Bj = Bq^ if coeffsupp(pj) — coeffsupp(qj), |coeffsupp(qj)| > 2, and one of coeff(qj, fc) 
and coeff(qj, — fc) is nonzero; otherwise Bj = Is. Note that Bo • • • Bj_i is paraunitary 
and satisfies 5(Bo---Bj_i) = {S9)*S6 with S9 = [1^^ , -1^2, z^^l^,, -z"^!^^]. By 
item (1) of Lemma [H the coefficient support of Bo • • • Bj_iBj is also contained in- 
side [—1,1]. By induction, the coefficient support of Bi • • • B^ is contained inside 
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[—1,1]. Moreover, Bi • • • takes the form in (|4.3p . Next, since Bi^_i; i;^ is con- 
structed recursively from pairs (qi,q2) of Qg Q(Bi • • • B^), by applying induction 
again and using item (2) of Lemma [1] we conclude that the coefficient support of 
B := (Bi • • • Bj.)B(-_fc is contained inside [—1, 1]. 



B 



Due to the Property (PI), (P2) of Bq and (P3), (P4) of B(q^^q,), Bi, . . . , B^ and 



(-fc,fc) 



reduce Q of the form in ^Ji) to Qi = Q(Bi • • • Br)B^ 



-k,k) 



QB of the form 



in (j2.8p with at least one of coefF(Qi,— fc) and coefF(Qi,fc) being 0. As constructed 
in section 2, Bq^ = Ig for the case that coeff(Qi,— fc) = coeff(Qi,fc) = 0, or Bq^ = 
diag([/iWi, /^g+sji; for the case coeff(Qi, fc) ^ 0, or Bq^ := <:\\a.g{Is^+s2,U3^s)E for 
the case that coeff(Qi, — fc) ^ 0. We next show that cocffsupp(BBQ^) C [—1, 1]. 

Let Q take the form in ()2.7p and Qi take the form in p.8|) with cocfF(Qi,fc) ^ 
0. Then Bq^ := diag(J7iWi, /s3+s^)£^ with C/i, Wi, and E being constructed as in 
section 2. Note that B takes the form in (14. 3p . Define 



[Gi, G2, G^jGqjFy, Fs] 



Gil G21 F'si F41 G51 Gai F71 Fsi 

G12 G22 F32 F42 G52 G62 -FV2 ^82 



By coeff(Qi, fc) = coeff(Q, fc - l)coeff(B, 1) + coeff(Q, fc)coefr(B, 0), we have 



Gi = G5A1 + G6A2 + F7B1 - F8B2 + GiDi + F3B1 + FiB2\ 
G2 = GsAa + G&A4. - F7C1 + F8C2 + G2D2 + F3C1 + F4G2; 

= Ft As + Fg^e + G1B3 + G2G3 + F3D3 F3; 

= F7A7 + FsAs + G1B4 + G2G4 + F4D4 — h, 



(4.5) 



where Gi,G2 are matrices defined in ()2.11|) . Then Ui = dia,g{U q^,Uq^) and Wi is 
defined as in (I2.12p . By the coefficient supports of B and Bq^ , we only need to check 
that coeff(B diag(C/iWi, /S3+S4), — 2) = 0. Let ^12, ^21, ^2 be diagonal matrices 



of size si X si, si x S2, S2 x si, S2 x S2, respectively, and satisfy diag(Vj< 
for j,£ ~ 1,2, where mi is the rank of Gi. Then 



[lmi,0] 



coeff(Bdiag(;7iWi, 1^3+34), -2) coeff(B, -1) • coeff(diag(J7iWi, 



S3 + Si)i 



-1) 



Ai 
-A2 







-A3 
A4 





B3 
-G3 

As 



-Bi 
Ci 

-A, 
Ag 



UqV21 UqV22 







Thus, we need to show AiU^^Vij — A^U^^V2j = and A2UQ^Vij — A4U^^V2i ~ 0, 
for j — 1,2, which is equivalent to showing that VjiU- Al — Vj2?7~ A^ = and 

VjiUt^Al - Vj2U*~^Al = for j = 1,2. Since GiC/g^ = [i?, 0] and G2U^^ = [i?,0], 
for some lower triangular matrix R of full rank mi, it is equivalent to proving that 
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GiAl - G2AI = and GiA* - = 0. By we have, 



= (G5A1 + G6^2 + FtBi - + GiDi + F3B1 + KiB2)At 

- {G5A3 + G6A4 - F7G1 + FsC2 + G2D2 + F3G1 + F4C2)A; 
+ (FtAs + FaAe + G1B3 + G2G3 + FiD-i)Bl 

- (FrAr + FgAg + G1B4, + G2C4 + FADi)Bl 
= Gr^AiAl ~ A^Al) + G^{A2A{ - ^4^;) 

+ Fj{BxAl + CxAl + A5B; - A7BI) 

+ Fs{-B2Al - G2AI + Ae,B; - AsBl) 

+ GiiDiAl + BsB; - B4BI) + G2{-D2A*3 + G3S3* - G4B4*) 

+ FaiBiAl - G1A3 + D3B3) + ^4(52^1 - G2^3 - D^Bl) = 0, 

where the last identity foUows from BB* = Is and coefF(QB, fc + 1) = 0. Similarly, 
G1A2 — G2AI = 0. The computation for showing coeffsupp(BBQ J C [—1,1] with 
Bq^ = diag(/s^+s2, L/3W3)i? is similar. Consequently, coeffsupp(BBQ J C [—1,1]. 
Therefore, item (1) holds. Item (2) is due to the property (|4.2I) of Uf and Ug- 

Note that SB = {Se)*Se with S0 = [1^,, z^^l^g, -z^^l^i]- And by the 
construction of Bq^, SBq^ — (50)* [Is'^ , — , z^^l^^ , — z~^lsj for some nonnegative 

integers s[, . . . , s'^ depending on the rank of Gi or G3 (see section 2). Consequently, 
item (3) holds. This also completes the proof of Algorithm [T] □ 
Now, we are ready to prove Theorems [1] and [H 

Proof of Theoremsl^and\^' The sufficiency part of Theorem[2]is obvious. We only 
need to show the necessary part. Suppose SP — {S9i)*S92. Let Q U^g^PU^gj and 
coeffsupp(Q) :— [fci,fc2]- Then SQ satisfies p.ip . By Lemma[2l the step of support 
reduction in Algorithm [1] produces a sequence of paraunitary matrices Ai,...,Aj 
with coefhcient support contained inside [—1,1] such that QAi • • • Aj = [/r, 0]. Due 
to item (1) of Lemma [2 J ^ f^^^r^l- Let Pj := A*, Po := U|g^ and P.7+1 := 
diag(U5ei,/s-r)- Then P^ := Pj+iPj • • • PiPq satisfies [/,.,0]Pe = P.^By item (3) of 
Lemma[21 (Pj+i, Pj) has mutually compatible symmetry for all ^ j ^ J. The claim 
that |coeffsupp([Pe]fcj)| ^ maxi^„^r |coeffsupp([P]„j)| for 1 ^ j,k ^ s follows from 
item (2) of Lemma [H Hence, all claims in Theorems [1] and [D have been verified. □ 
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